Sharp phase transitions in a small frustrated network of trapped ion spins 
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Sharp quantum phase transitions typically require a large system with many particles. Here we 
show that for a frustrated fully-connected Ising spin network represented by trapped atomic ions, 
the competition between different spin orders leads to rich phase transitions whose sharpness scales 
exponentially with the number of spins. This unusual finite-size scaling behavior opens up the 
possibility of observing sharp quantum phase transitions in a system of just a few trapped ion spins. 
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Quantum simulators are motivated by the promise of 
gaining insight into many-body quantum systems such 
as high- Tc superconductors or complex arrangements of 
interacting spins. Cold atomic systems form a promising 
platform for quantum simulation, as the interactions be- 
tween particles can be under great control. A good exam- 
ple is a collection of trapped and laser-cooled atomic ions, 
each representing an effective spin that can be made to 
interact with all the others by modulating the Coulomb 
interaction between ions. By applying spin-dependent 
optical dipole forces it has been shown that a crystal of 
ions provides an ideal platform to simulate intractable 
interacting spin models 0, Q ■ Following this proposal, 
recent experiments have simulated quantum magnetism 
with a few ions For three or more ions, the long- 

range coupling between the spins can provide frustrated 
interaction patterns or competition between various spin 
orders and interesting quantum phase diagrams p| [|| . 
However, the observation of a quantum phase transition 
typically requires a large system with many particles, as 
the width (sharpness) of a quantum phase transition usu- 
ally scales with 1/N, the inverse of the number of par- 
ticles [7|. With such slow finite-size scaling laws, small 
networks of trapped ions realized in current experiments 
(N < 20) are not expected to exhibit sharp transitions 
between distinct quantum phases. 

In this paper, we show the surprising result that sharp 
phase transitions can indeed be observed with just a few 
atomic ions. This is due to unusual finite size scaling 
laws in this frustrated spin network, where the sharpness 
of some phase transitions scales exponentially instead of 
linearly with 1/N. By controlling a single experimen- 
tal parameter that determines the pattern of spin-spin 
couplings between the ions, we show that the expected 
ground state emerges from a delicate compromise be- 
tween the couplings. Frustration in the spin network 
leads to a variety of spin orders, with the number of dis- 
tinct phases increasing rapidly with the number of ions. 
We construct the complete phase diagram for small spin 
networks realizable with the current technology. The 
sharp phase transition is characterized in detail with an 
explanation of its unusual finite-size scaling behavior. 

We consider a small crystal of ions confined in a 



one-dimensional harmonic trap. The spin states of the 
ions are represented by two internal states, referred as 
|t) and |4_), and the effective spin-spin interaction be- 
tween the ions is induced with off-resonant bichromatic 
laser beams 0, [^-Q. The ion-laser coupling Hamilto- 
nian, written in the rotating frame, has the form H = 
J2 n [fj£lcos(5kx n + yA)a z n + Hi where fi is a Ra- 

man Rabi frequency, 8k is the wave vector difference be- 
tween the two Raman beams (which is assumed to be 
along the radial direction x), [i is the beatnote or detun- 
ing between the two laser beams, and are Pauli ma- 
trices describing the spin of the nth ion, and B is an effec- 
tive magnetic field induced by radiation that coherently 
flips the spins. In the rotating frame, the radial coordi- 
nate x n is expanded in terms of the transverse phonon 
modes a k as x n = J2 k &,Vft/(2mw fc )(a*e^ f +a t e- KJ '= i ), 
where m is the atomic mass, Uk is the eigenfrequency of 
the fcth normal mode of the ion crystal, and b„ is the 
eigenmode transformation matrix. We use transverse 
phonon modes because they can more easily be scaled 
up to large systems [a, |9( . Under the Lamb-Dicke crite- 
rion r] ni k = b^ l 5ky/h/(2muk) <C 1, the Hamiltonian H 

is simplified to H = — hfl J2 n k k sm (/•**) a n( a l elUkt + 
a fc e-*"**) + fl£«<7* 

If we assume that the laser detuning u is not resonant 
with any phonon mode with the condition \u>k — n\ 3> 
rjn^Q, satisfied for all n modes k, the probability of ex- 
citing any phonon mode \Qr) Tlt kb^ l /2 (u>k — /i)| 2 is negligi- 
ble. We can therefore adiabatically eliminate the phonon 
modes and arrive at the following effective spin-spin cou- 
pling Hamiltonian 0, [13] 



Hf, — > J, 



where the coefficients 



J m,n. — 



{httsky 

2m 



E 



b k b k 



CO 



2 ' 



(1) 



(2) 



This Ising Hamiltonian is a pillar of many-body 
physics, and its properties have been exhaustively stud- 
ied under various conditions Q- For instance, the ground 
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Figure 1: (Color online) Illustration of spin orders in two frus- 
trated Ising networks with competing long-range interaction 
for N = 7 ions. Left panel: Ferromagnetic (FM) order with 
detuning fi = 5.1; Right panel: Kink order with fi = 5.3. The 
thickness of the edge in the graph represents the strength of 
each coupling. Positive (AFM) spin couplings are indicated 
in red, while negative (FM) couplings are indicated in black. 



state of the Ising Hamiltonian is well understood when 
the coupling coefficients J mn are uniform, or nonzero 
only for nearest neighbors. However, here we have an ex- 
tended Ising network where the coupling coefficients J mn 
are inhomogeneous (both in magnitude and sign) and ex- 
tend over long range [ll[ . The strong competition among 
these interaction terms (even with B — 0) will generally 
lead to highly frustrated ground states where individual 
bonds are compromised in order to reach a global en- 
ergy minimum. For arbitrary coupling coefficients J mn , 
the Hamiltonian (1) generally belongs to the complexity 
class of NP-complete problems [12j , meaning that calcu- 
lating attributes such as the system ground state becomes 
intractable when the system size is scaled up. 

We consider the case where the coupling coefficients 
J mn are controlled by a single experimental parameter, 
the laser detuning /j, To determine J mn from any 

detuning fi with the formula (2), we first need the normal 
mode transformation matrix b\. This is obtained by find- 
ing the equilibrium positions for a given number of ions in 
a harmonic trap and then diagonalizing the Coulomb in- 
teraction Hamiltonian expanded about the ions' equilib- 
rium positions. With a single control parameter fi, we are 
not able to program arbitrary coupling coefficients J m „. 
However, the interaction pattern is sufficiently complex 
to allow frustrated ground state configurations and rich 
phase transitions. To illustrate this we show in Fig. 1 a 
coupling pattern for N = 7 ions and its associated ground 
state spin configuration for B — 0. The coupling pattern 
is represented by a graph where the color and the thick- 
ness of each edge represents respectively the sign (ferro- 
magnetic or antiferromagnetic) and the magnitude of the 
coupling. In Fig. 1(a), we find a ferromagnetically or- 
dered ground state with all the spins pointing to the same 
direction. However, in this ferromagnetic state, some of 
the bonds, such as the strong antiferromagnetic bond be- 
tween the ions 1 and 7, are compromised, and due to this 
frustration, the ground-state spin configuration is very 
sensitive to the strength of the coupling. If we adjust the 
detuning fi by a small fraction of the trap frequency, the 
ferromagnetic bonds of the ion pairs (1,5) and (3,7) are 



slightly weakened (see Fig. 1(b)) and the antiferromag- 
netic bond (1,7) dominates and flips the spin direction 
of the entire left (or right) half of the ion crystal. This is 
a phase transition from ferromagnetic order to a "kink" 
order, with a kink in the spin direction between the 4th 
and 5th ions counting from either the left or the right 
side. 

To show the rich phase diagram for this system, in Fig. 
2(a) we list all different spin phases at B — for a small 
Ising network with 3, 5, 7 and 9 ions. For an odd number 
of ions, the phase diagram is more interesting and fea- 
tures a larger variety of spin orders, because the left-right 
reflection symmetry in a linear ion crystal can be sponta- 
neously broken. In Fig. 2, for convenience, we denote the 
laser detuning fi with a re-scaled dimensionless parameter 
ft by labeling the phonon-mode eigen-frequencies (from 
lowest to highest) with integers (from 1 to N for an iV-ion 
crystal) and using a linear scaling between two integers to 
denote detuning located between the corresponding two 
modes. For instance, in this notation // = 2.75 means the 
physical detuning fi = uji + 0.75(w3 — u)%). Each phase 
is characterized by a spin order (denoted with a binary 
string where and 1 correspond to t and 4- spin respec- 
tively) which gives one of the ground state spin configu- 
rations. The Ising Hamiltonian (1) features a reflection 
symmetry and an intrinsic Z2 symmetry with respect to 
a global spin flip. The spin order breaks the Ising symme- 
try, so each phase is at least two-fold degenerate. If the 
spin order also breaks the reflection symmetry, the corre- 
sponding ground state is 4-fold degenerate. For instance, 
for the phase denoted by the spin order 01001, the four 
degenerate ground states are ||t44t), |T4-Tt4-)> |t44-t4-)> 
and |4-tt4t)- When the re-scaled detuning f~i crosses an 
integer (a phonon mode), the spin order changes as ex- 
pected, but this is not a conventional phase transition as 
the parameters J mn change discontinuously in the Hamil- 
tonian (1). However, when fi varies within two adjacent 
integers, all the parameters J mn are analytic functions 
of jl, yet the spin order can still change abruptly, sig- 
naling a phase transition. The frequency of this type of 
inter-mode phase transition increases rapidly with the ion 
number: there is one such transition for a three-ion chain 
and 12 such transitions in a nine-ion crystal. Another no- 
table feature from Fig. 2(a) is that there is typically no 
phase transition when fi varies from an even mode (2A;) 
to an odd mode (2k + 1, k = 0, 1,- • • , (N - 1) /2). In 
such regions, the spin order has a reflection symmetry. 
This suggests that a spin order with reflection symmetry 
may be more stable in energy and does not easily yield 
to other spin configurations. This observation is consis- 
tent with the fact that for an even number of ions, there 
are much fewer inter-mode phase transitions, as the spin 
order in these cases has a reflection symmetry. 

As we add a transverse B field to the Hamiltonian, 
the spins will gradually become polarized along the in- 
direction along B. In Fig. 2(b), we plot the aver- 
age polarization (^2 n &n)/N as a function of the field 
B (in the unit of the average J mn defined by J = 
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Figure 2: (Color online) (a) Ground-state phase diagram at B — characterized by the corresponding spin orders for N = 
3,5,7,9 ions. The transition point is denoted by the re-scaled detuning ft. (b) Average polarization (y2 n c?i^/N for N — 7 
ions at finite fields B. The right figure is a closeup near the critical point. 



^E m ^n\Jmn\ 2 /[N(N-l)}) for TV = 7 ions in a small 

region of the detuning ft. We find that the system is eas- 
ily polarized if it lies at the critical point between two 
different spin orders given by the Ising couplings. But 
near the center of a spin phase, the spin order is more ro- 
bust and can persist under a finite B, eventually yielding 
to the polarized phase as B increases through the Ising- 
type transition (which becomes a broad crossover for this 
finite system). 

With B = 0, the transition between different spin or- 
ders is sharp as it is characterized by a level crossing 
for the ground state of the Hamiltonian (1). When we 
turn on a finite B field, the system shows only avoided 
level crossings in its ground state, and we expect that 
the sharp phase transition at B = to be replaced by a 
broad crossover for this small system, similar to the Ising- 
type of transition discussed above. Interestingly, this is 
not always the case. Even at a finite B, some transitions 
remain very sharp in this small system, and the sharp- 
ness increases exponentially with the number of particles. 
This provides an unusual finite-size scaling behavior, dis- 
tinctively different than the polynomial finite-size scaling 
of the transition width that one typically sees in conven- 
tional phase transitions [7]. (For instance, the conven- 
tional Ising transition displays a linear scaling law (7|.) 
To characterize this unusual finite-size scaling behavior, 



we look at a particular example: for N ions, there is a 
unique spin phase transition in the region with detuning 
N — 2 < ft < N — 1 when iV is odd. A schematic phase 
diagram for this region is shown in Fig. 3(a). At B = 0, 
we have a ferromagnetic phase on the left side which is 
doubly degenerated and a kink phase on the right side 
which is 4-fold degenerate. At finite B, these two spin 
orders remain robust in a range of B, before eventually 
yielding to a polarized phase for a large B field through 
a crossover. The transition between the ferromagnetic 
and the kink phases remain sharp for a finite B, and 
this sharpness increases rapidly with the ion number, as 
one can see from Fig. 3(b) and 3(c). With a moderate 
increase of the ion number from 5 to 9, the transition 
between the ferromagnetic and the kink phase becomes 
almost infinitely sharp. It is also interesting to note that 
the transition line between these two phases has a slope 
with the B axis, so one can cross this phase transition by 
tuning either the detuning ft or the field B. 

We have characterized the transition width between 
the ferromagnetic phase and the kink phase, and find this 
transition width indeed shrinks exponentially with the 
number of ions. To see this, we show the four lowest lev- 
els across this phase transition when we tune either ft or 
B (see Fig. 4a) . While the lowest state is a smooth func- 
tion of ft at a finite B, the second and third states feature 
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Figure 3: (Color online) (a) The schematic phase diagrams for odd numbers of ions in the region with detuning TV— 2 < p, < TV— 1. 
For a small number of ions, the solid line represents a sharp transition, whereas the dashed lines represent a continuous crossover 
to the polarized state. (b,c) The calculated theoretical phase diagrams for (b) TV = 5 and (c) TV = 9 ions. Color shows the 
order parameter defined by Pfm — Pk, where Pfm and Pk are the projection probabilities of the ground state of the system to 
the Hilbert subspace with the ferromagnetic and the kink orders, respectively (the basis-vectors of the corresponding subspaces 
are indicated in Fig. (a)). 




Figure 4: (a) The structure of the lowest four energy levels 
around the sharp ferromagnetic-kink phase transition. The 
two lowest states have a ferromagnetic (kink) order on the 
left (right) side. The transition width W is proportional to 
the energy gap AE, and AE is fitted by an exponential, (b) 
The exponent a is a function of the ion number TV, where the 
solid line is a linear fit a ~ (TV — 1) /2. 

a level crossing. From Fig. 4a, the transition width W 
is proportional to the energy gap AE. We calculate the 
energy gap AE as a function of the ion number TV, and 

find it is well fit by AE cx (b/NJ) " for B < TV J, where 

the exponent a has a roughly linear dependence with TV 
as shown in Fig. 4b and can be fit by a ~ (TV — 1) /2. 
This formula shows that the energy gap AE, and thus 
also the transition width W, shrinks exponentially with 



the ion number, which explains why the transition is so 
sharp for a small system. The exponential shrinking of 
the energy gap AE can be intuitively understood as fol- 
lows: when B <§C TV J we can treat the term B^2 n a^ l 
as a perturbation in the Hamiltonian (f). For each ap- 
plication of B^2 n (Tni we can only flip the direction of 
one spin. As the ferromagnetic state and the kink state 
have (TV — 1) /2 spins taking opposite directions, the two 
states need to be connected through (TV — 1) /2th order 
perturbation, and thus the energy gap is proportional to 

(B/NJ) (N - 1)/2 . 

In summary, we have shown that laser induced mag- 
netic coupling between trapped ions realizes a frustrated 
Ising spin network with competing long range interac- 
tions, giving rise to rich phase diagrams for the ground 
state. Some of the phase transitions in this system are 
characterized by a unusual finite size scaling, where the 
transition width scales down exponentially with the num- 
ber of ions. This exponential finite size scaling lads to 
sharp phase transitions for a small system even with just 
a few ions, as one can realize now in the lab. 
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